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The crustal fraction of moment of inertia in neutron stars is calculated using /3-equilibrated nuclear 
matter obtained from density dependent M3Y effective interaction. The transition density, pressure 
and proton fraction at the inner edge separating the liquid core from the solid crust of the neutron 
stars determined from the thermodynamic stability conditions are found to be pt = 0.0938 fm~®, 

Pt = 0.5006 MeV fm“® and Xp^t) = 0.0308, respectively. The crustal fraction of the moment of 
inertia can be extracted from studying pulsar glitches and is most sensitive to the pressure as well 
as density at the transition from the crust to the core. These results for pressure and density at 
core-crust transition together with the observed minimum crustal fraction of the total moment of 
inertia provide a new limit for the radius of the Vela pulsar: A > 4.10 + 3.SGM/Mq kms. 
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I. INTRODUCTION 

Pulsar glitches, which are discontinuities in the spin- 
down of pulsars, involve sudden transfer of angular mo¬ 
mentum from an isolated component (consisting of super¬ 
fluid neutrons in crust) to the entire star through vortex 
unpinning. The sudden jumps in rotational frequencies 
w which may be as large as ^ ^ 10“® — 10“® have 
been observed for many pulsars. The frequency of ob¬ 
served glitches is statistically consistent with the hypoth¬ 
esis that all radio pulsars experience glitches [l| . Glitches 
are thought to originate from interactions between the 
rigid neutron star crust, typically somewhat more than 
a kilometer thick, and rotational vortices in a neutron 
superfluid. The inner kilometer of crust consists of a 
crystal lattice of nuclei immersed in a neutron superfluid 
Q. Because the pulsar is spinning, the neutron super¬ 
fluid (both within the inner crust and deeper inside the 
star) is threaded with a regular array of rotational vor¬ 
tices. As the spin of the pulsar gradually slows, these 
vortices must gradually move outwards since the rota¬ 
tional frequency of a superfluid is proportional to the 
density of vortices. Deep within the star, the vortices 
are free to move outwards. In the crust, however, the 
vortices are pinned by their interaction with the nuclear 
lattice. Various theoretical models differ in impor¬ 
tant respects as to how the stress associated with pinned 
vortices is released in a glitch: for example, the vortices 
may break and rearrange the crust, or a cluster of vor¬ 
tices may suddenly overcome the pinning force and move 
macroscopically outward, with the sudden decrease in the 
angular momentum of the superfluid within the crust re¬ 
sulting in a sudden increase in angular momentum of the 
rigid crust itself and hence a glitch. All the models agree 
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that the fundamental requirements are the presence of 
rotational vortices in a superfluid and the presence of a 
rigid structure which impedes the motion of vortices and 
which encompasses enough of the volume of the pulsar to 
contribute significantly to the total moment of inertia. 

In the present work, the equation of state (EoS) used 
is obtained from the density dependent M3Y effective 
nucleon-nucleon interaction (DDM3Y) for which the in¬ 
compressibility Kao for the symmetric nuclear matter 
(SNM), nuclear symmetry energy EgymiPo) at satura¬ 
tion density po, the isospin dependent part Kp of the iso- 
baric incompressibility and the slope L are in excellent 
agreement with the constraints recently extracted from 
measured isotopic dependence of the giant monopole res¬ 
onances in even-A Sn isotopes, from the neutron skin 
thickness of nuclei, and from analyses of experimental 
data on isospin diffusion and isotopic scaling in interme¬ 
diate energy heavy-ion collisions [ 1 , The core-crust 
transition in neutron stars is determined 0 by analyz¬ 
ing the stability of the /3-equilibrated dense nuclear mat¬ 
ter with resp ect to the thermodynamic stability condi¬ 
tions |ll| - [l^ . The mass-radius relation for neutron stars 
is obtained by sol ving the Tolman-Oppenheimer-Volkoff 
Equation (TOV) [Ig, [I4I and then the crustal fraction 
of moment of inertia is determined using pressure and 
density at core-crust transition. As the angular momen¬ 
tum requirements of glitches in Vela pulsar indicate that 
1.4% of the star’s moment of inertia drives these events, 
the allowed region for masses and radii for Vela pulsar is 
determined from the condition that the crustal fraction 
of moment of inertia ^ > 0.014 which sets a new limit 
for its radius. 


II. CORE-CRUST TRANSITION IN 
/3-EQUILIBRATED NEUTRON STAR MATTER 

The nuclear matter EoS is calculated using the 
isoscalar and the isovector [l^, [l^ components of M3Y 
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interaction along with density dependence. The density 
dependence of this DDM3Y effective interaction is com¬ 
pletely determined from nuclear matter calculations. The 
equilibrium density of the nuclear matter is determined 
by minimizing the energy per baryon. The energy varia¬ 
tion of the zero range potential is treated accurately by 
allowing it to vary freely with the kinetic energy part 
gfem qI? energy per baryon e over the entire range of e. 
This is not only more plausible, but also yields excellent 
result for the incompressibility K^o of the SNM which 
does not suffer from the superluminosity problem [^ . 
In a Fermi gas model of interacting neutrons and pro¬ 
tons, with isospin asymmetry X = , p = Pn + Pp, 

where pn, Pp and p are the neutron, proton and baryonic 
number densities respectively, the energy per baryon for 
isospin asymmetric nuclear matter can be derived as [2^ 
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where nib is the baryonic rest mass, kp={1.57T^p)^ 
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and Jvoi represent the volume integrals of the isoscalar 
and the isovector parts of the M3Y interaction. The 
isoscalar and the isovector components of 

M3Y interaction potential are given by 


ing hciSn^pXp)^^^ = —= +2-^. The thermody¬ 
namical method requires the system to obey the intrinsic 
stability condition Vthermai > 0 which is given by 
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and obviously, it goes to zero at the inner edge separating 
the liquid core from the solid crust since it corresponds to 
a phase transition from the homogeneous matter at high 
densities to the inhomogeneous matter at low densities. 
The core-crust transition density pt, pressure Pt and pro¬ 
ton fraction Xp(^t) of the neutron stars are obtained [loj| 
by setting Vthermai = 0 which goes to negative with de¬ 
creasing density. 


III. CRUSTAL FRACTION OF MOMENT OF 
INERTIA IN NEUTRON STARS 

The crustal fraction of the moment of inertia ^ can 
be expressed as a function of M (gravitational mass of 
the star) and R (radius of the star) with the only depen¬ 
dence on the equation of state arising from the values of 
transition density pt and pressure P*. Actually, the ma¬ 
jor dependence is on the value of Pt, since pt enters oidy 
as a correction in the following approximate formula 
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where s represents the relative distance between two in¬ 
teracting baryons, Jqo = —276 MeV.fm^, Jqi = -1-228 
MeV.fm^ and the energy dependence parameter a = 
0.005 MeV“^. The strengths of the Yukawas were ex¬ 
tracted by fitting its matrix elements in an oscillator basis 
to those elements of G-matrix obtained with Reid-Elliott 
soft core NN interaction and the ranges were selected to 
ensure OPEP tails in the relevant channels as well as a 
short-range part which simulates the cr-exchange process 
[2l|. The density dependence is employed to account for 
the Pauli blocking effects and the higher order exchange 
effects [ 2 ^. Thus the DDM3Y effective NN interaction 
is given by uoi(s, p, e^*") = Wp) where the 

density dependence g{p) = (7(1 — /3p") 0 with C and 
/3 being the constants of density dependence. 

The /3-equilibrated nuclear matter EoS is obtained by 
evaluating the asymmetric nuclear matter EoS at the 
isospin asymmetry X = 1 — 2xp determined from the /3- 
equilibrium proton fraction Xp [= ^], obtained by solv- 
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where The crustal fraction of the moment of 

inertia is particularly interesting as it can be inferred 
from observations of pulsar glitches, the occasional dis¬ 
ruptions of the otherwise extremely re gula r pulsations 
from magnetized, rotating neutron stars |23j . Link et al. 
Q showed that glitches represent a self-regulating insta¬ 
bility for which the star prepares over a waiting time. 
The angular momentum requirements of glitches in the 
Vela pulsar indicate that more than 0.014 of the moment 
of inertia drives these events. So, if glitches originate in 
the liquid of the inner crust, this means that ^ > 1.4%. 


IV. TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATION AND MASS-RADIUS RELATION 

In astrophysics, the Tolman-Oppenheimer-Volkoff 
(TOV) equation ul, [13 constrains the structure of a 
spherically symmetric body of isotropic material which is 
in static gravitational equilibrium, as modelled by gen¬ 
eral relativity and is given by 



























FIG. 1: Variation of mass with central density for slowly ro¬ 
tating neutron stars for the present nuclear EoS. 
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which can be easily solved numerically using Runge- 
Kutta method for masses and radii. The quantities e{r) 
and P{r) are the energy density and pressure at a radial 
distance r from the centre, and are given by the equa¬ 
tion of state. The mass of the star contained within a 
distance r is given by m{r). The size of the star is de¬ 
termined by the boundary condition P{r) = 0 and the 
total mass M of the star integrated up to the surface R is 
given by M = m{R) [l^. The single integration constant 
needed to solve the TOV equation is Pc, the pressure at 
the center of the star calculated at a given central density 
Pc- The masses of slowly rotating neutron stars are very 
close [2^ - 1^ to those obtained by solving TOV equation. 

The moment of inertia of neutron stars is calculated by 
assuming the star to be rotating slowly with a uniform 
angular velocity [1^. The angular velocity uj{r) of a 
point in the star measured with respect to the angular 
velocity of the local inertial frame is determined by the 
equation 
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dr 
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FIG. 2: The mass-radius relation of slowly rotating neutron 
stars for the present nuclear EoS. For the Vela pulsar, the 
constraint of -^ > 1.4% implies that allowed masses and 
radii lie to the right of the line defined by -^ = 0.014 (for 
pt = 0.0938 fm"®, Ft = 0.5006 MeV fm"®). 


The function (j){r) is constrained by the condition 


p{r) = constant = /i(P) 


r 2GM 
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where the chemical potential p{r) is defined as 
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Using these relations, Eq.(6) can be solved subject to the 
boundary conditions that w(r) is regular as r —>■ 0 and 
uj{r) —>■ U as r —^ oo. Then moment of inertia of the star 
can be calculated using the definition I = J/U, where 
the total angular momentum J is given as 
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V. THEORETICAL CALCULATIONS 

The calculations are performed using the values of the 
saturation density po=0.1533 fm“^ and the satu¬ 
ration energy per baryon eo=-15.26 MeV [1^ for the 
SNM obtained from the co-efficient of the volume term of 
Bethe-Weizsacker mass formula which is evaluated by fit¬ 
ting the recent experimental and estimated atomic mass 
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TABLE I: Results of present calculations for u=| of symmet¬ 
ric nuclear matter incompressibility Koo, nnclear symmetry 
energy at saturation density Esym{po), the slope L and isospin 
dependent part Kt of the isobaric incompressibility (all in 
MeV) are tabulated along with the density, pressure and 
proton fraction at the core-crust transition for /?-equilibrated 
neutron star matter and corresponding Vela pulsar constraint. 



Esymi^po') 

L 

K-r 

274.7 ± 7.4 

30.71 ±0.26 

45.11 ±0.02 

-408.97 ±3.01 

? 

1 

CO 

Pt(MeVfm-3) 

Xp(t) 

Vela pulsar R(km) 

0.0938 

0.5006 

0.0308 

A > 4.10 ± 3.36M/Mq 


excesses from Audi-Wapstra-Thibault atomic mass table 
by minimizing the mean square deviation in corp orat- 
ing correction for the electronic binding energy |32j ]. In a 
similar recent work, including surface symmetry energy 
term, Wigner term, shell correction and proton form fac¬ 
tor correction to Coulomb energy also, turns out to 
be 15.4496 MeV and 14.8497 MeV when and 
terms are also included [s^. Using the usual value of 
a = 0.005 MeV“^ for the parameter of energy depen¬ 
dence of the zero range potential and n = |, the values 
obtained for the constants of density dependence C and 
P and the SNM incompressibility K^o are 2.2497, 1.5934 
fm^ and 274.7 MeV, respectively. The saturation en¬ 
ergy per baryon is the volume energy coefficient and the 
value of -15.26±0.52 MeV covers, more or less, the entire 
range of values obtained for av for which now the val¬ 
ues of C'=2.2497±0.0420, ^=1.5934±0.0085 fm^ and the 
SNM incompressibility Aroo=274.7±7.4 MeV [^ . 

The stability of the /3-equilibrated dense matter in 
neutron stars is investigated and the location of the in¬ 
ner edge of their crusts and core-crust transition density 
and pressure are determined using the DDM3Y effective 
nucleon-nucleon interaction. The results for the transi¬ 
tion density, pressure and proton fraction at the inner 
edge separating the liquid core from the solid crust of 
neutron stars are calculated and presented in Table-I for 
n = |. The symmetric nuclear matter incompressibil¬ 
ity Krxn nuclear symmetry energy at saturation density 
Esymipo), the slope L and isospin dependent part K-r 
of the isobaric incompressibility are also tabulated since 
these are all in excellent agreement with the constraints 
recently extracted from measured isotopic dependence of 
the giant monopole resonances in even-A Sn isotopes, 
from the neutron skin thickness of nuclei, and from anal¬ 
yses of experimental data on isospin diffusion and isotopic 
scaling in intermediate energy heavy-ion collisions. 

The calculations for masses and radii are performed us¬ 
ing the EoS covering the crustal region of a compact star 
which are Feynman-Metropolis-Teller (FMT) [s^, Baym- 
Pethick-Sutherland (BPS) [s^ and Baym-Bethe-Pethick 
(BBP) upto number density of 0.0582 fm“^ and /3- 
equilibrated neutron star matter beyond. The values of 
/ obtained by solving Eq.(6) subject to the boundary 
conditions stated earlier are listed in Table-II along with 


TABLE II: Radii, masses, total & crustal fraction of moment 
of inertia and crustal thickness as functions of central density 
Pc- 


Pc 

R 

M 

I 

AI 

T 

AR 

fm"® 

km 

Mq 

MQknP 

fraction 

km 

2.00 

8.6349 

1.8277 

70.88 

0.0055 

0.2462 

1.90 

8.7598 

1.8467 

73.83 

0.0057 

0.2523 

1.80 

8.8957 

1.8651 

77.00 

0.0060 

0.2598 

1.70 

9.0444 

1.8824 

80.38 

0.0063 

0.2696 

1.60 

9.2052 

1.8980 

83.97 

0.0067 

0.2806 

1.50 

9.3810 

1.9109 

87.70 

0.0072 

0.2951 

1.40 

9.5710 

1.9197 

91.52 

0.0079 

0.3121 

1.39 

9.5911 

1.9203 

91.91 

0.0080 

0.3144 

1.38 

9.6109 

1.9208 

92.29 

0.0080 

0.3161 

1.37 

9.6314 

1.9213 

92.67 

0.0081 

0.3185 

1.36 

9.6514 

1.9217 

93.05 

0.0082 

0.3203 

1.35 

9.6718 

1.9220 

93.43 

0.0083 

0.3222 

1.34 

9.6928 

1.9223 

93.81 

0.0084 

0.3248 

1.33 

9.7141 

1.9225 

94.18 

0.0085 

0.3275 

1.32 

9.7349 

1.9226 

94.55 

0.0085 

0.3296 

1.31 

9.7559 

1.9227 

94.93 

0.0086 

0.3318 

1.30 

9.7770 

1.9226 

95.30 

0.0087 

0.3340 

1.20 

9.9995 

1.9173 

98.85 

0.0098 

0.3620 

1.10 

10.2371 

1.9004 

101.88 

0.0112 

0.3970 

1.00 

10.4902 

1.8675 

103.94 

0.0132 

0.4441 

0.90 

10.7544 

1.8127 

104.42 

0.0158 

0.5066 

0.80 

11.0239 

1.7285 

102.47 

0.0197 

0.5929 

0.70 

11.2865 

1.6064 

97.04 

0.0255 

0.7148 

0.60 

11.5245 

1.4369 

87.06 

0.0344 

0.8952 

0.59 

11.5456 

1.4170 

85.78 

0.0356 

0.9175 

0.58 

11.5666 

1.3965 

84.44 

0.0368 

0.9411 

0.57 

11.5874 

1.3753 

83.04 

0.0381 

0.9663 

0.56 

11.6073 

1.3536 

81.58 

0.0394 

0.9924 

0.55 

11.6262 

1.3313 

80.07 

0.0408 

1.0193 

0.50 

11.7135 

1.2104 

71.65 

0.0492 

1.1792 

0.45 

11.7830 

1.0734 

61.88 

0.0602 

1.3897 

0.40 

11.8388 

0.9206 

51.00 

0.0752 

1.6801 

0.30 

12.0129 

0.5808 

28.54 

0.1249 

2.7618 

0.25 

12.3703 

0.4103 

19.24 

0.1686 

3.9149 

0.24 

12.5113 

0.3779 

17.73 

0.1805 

4.2542 

0.23 

12.6944 

0.3464 

16.35 

0.1942 

4.6511 

0.22 

12.9314 

0.3159 

15.14 

0.2103 

5.1189 

0.21 

13.2434 

0.2867 

14.12 

0.2296 

5.6802 

0.20 

13.6576 

0.2587 

13.31 

0.2537 

6.3643 

0.19 

14.2131 

0.2323 

12.74 

0.2847 

7.2125 

0.18 

14.9725 

0.2075 

12.47 

0.3265 

8.2904 

0.17 

16.0398 

0.1845 

12.59 

0.3863 

9.7057 

0.16 

17.5771 

0.1634 

13.25 

0.4767 

11.6254 

0.15 

19.8913 

0.1445 

14.77 

0.6254 

14.3634 

0.14 

23.5740 

0.1278 

17.88 

0.8972 

18.5215 
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TABLE III: Variations of the core-crust transition density, pressnre and proton fraction for /3-equilibrated neutron star matter, 
symmetric nnclear matter incompressibility Kcx> and isospin dependent part AV of isobaric incompressibility with parameter n. 


n 

pt 

Pt 

Xp(t) 


Kr 

Maximum mass 

Radius 

Crustal Thickness 

Vela Pulsar Radius Constraint 



MeVfm-® 


MeV 

MeV 

Mq 

km 

km 

km 

Expt. 

values 

- 


250-270 

-370±120 

1.97±0.04 




1/6 

0.0797 

0.4134 

0.0288 

182.13 

-293.42 

1.4336 

8.5671 

0.4009 

R > 4.48 4- 3.37M/Mo 

1/3 

0.0855 

0.4520 

0.0296 

212.98 

-332.16 

1.6002 

8.9572 

0.3743 

R > 4.31 -h 3.36M/M0 

1/2 

0.0901 

0.4801 

0.0303 

243.84 

-370.65 

1.7634 

9.3561 

0.3515 

R > 4.19-h3.36M/M0 

2/3 

0.0938 

0.5006 

0.0308 

274.69 

-408.97 

1.9227 

9.7559 

0.3318 

R > 4.1O-h3.36M/M0 

1 

0.0995 

0.5264 

0.0316 

336.40 

-485.28 

2.2335 

10.6408 

0.3088 

R > 3.99 -h 3.36M/M0 


masses M, radii R and crustal thickness AR of neutron 
stars. Once masses and radii are determined, ^ are ob¬ 
tained from Eq.(4) and listed in Table-II. In Fig.-l, vari¬ 
ation of mass with central density is plotted for slowly 
rotating neutron stars for the present nuclear EoS. In 
Fig.-2, the mass-radius relation of slowly rotating neu¬ 
tron stars is shown. Using Eq.(4) again the mass-radius 
relation is obtained for fixed values of pt and Pj. 
This is then plotted in the same figure for ^ equal to 
0.014. For Vela pulsar, the constraint ^ > 1.4% implies 
that allowed mass-radius lie to the right of the line de¬ 
fined by ^ = 0.014 (for pt = 0.0938 fm-^, P* = 0.5006 
MeV fm“^). This condition is given by the inequality 
R > 4.10 -b 3.36M/Mq kms. 

In Table-Ill, the variations with parameter n of the 
core-crust transition density, pressure and proton frac¬ 
tion for /3-equilibrated neutron star matter, symmetric 
nuclear matter incompressibility iToo, isospin dependent 
part Kr of isobaric incompressibility, neutron star’s max¬ 
imum mass with corresponding radius and crustal thick¬ 
ness are listed along with corresponding Vela pulsar con¬ 
straints. It is important to mention here that recent ob¬ 
servations of the binary millisecond pulsar JI6I4-2230 
by P. B. Demorest et al. suggest that the masses lie 
within 1.97±0.04 Mq where Mq is the solar mass. 


VI. RESULTS AND DISCUSSION 

Recently, it is conjectured that the glitches observed in 
the Vela pulsar require an additional reservoir of angular 
momentum and the crust may not be enough to explain 
the phenomenon [^ . Large pulsar frequency glitches 
can be interpreted as sudden transfers of angular mo¬ 
mentum between the neutron superfluid permeating the 
inner crust and the rest of the star. In spite of the ab¬ 
sence of viscous drag, the neutron superfluid is strongly 
coupled to the crust due to non-dissipative entrainment 
effects. It is often argued that these effects may limit 
the maximum amount of angular momentum that can 
possibly be transferred during glitches [s^ . We find that 
the present EoS can accommodate large crustal moments 
of inertia and that large enough transition pressures can 


be generated to explain the large Vela glitches without 
invoking an additional angular-momentum reservoir be¬ 
yond that confined to the solid crust. Our results suggest 
that the crust may be enough [i^ which can be substan¬ 
tiated from Table-II that ^ > 0.014 for pulsars with 
masses 1.8 Mq or less. 

The results listed in Table-Ill suggest that SNM incom¬ 
pressibility do have some effect in determining the crustal 
fraction of moment of inertia and on the Vela Pulsar Ra¬ 
dius Constraint like some other recent studies [4l|. But 
the incompressibility values of about 15 MeV window 
around 270 MeV corresponding to n = | is experimen¬ 
tally supported. The present status of experimental de¬ 
termination of the SNM incompressibility from the com¬ 
pression modes ISGMR and isoscalar giant dipole reso¬ 
nance (ISGDR) of nuclei infers [l^ that due to violations 
of self consistency in HF-RPA calculations of the strength 
functions of giant resonances result in shifts in the calcu¬ 
lated values of the centroid energies which may be larger 
in magnitude than the current experimental uncertain¬ 
ties. In fact, the prediction of Koo lying in the range 
of 210-220 MeV were due to the use of a not fully self- 
consistent Skyrme calculations [d^. Correcting for this 
drawback, Skyrme parmetrizations of SL y4 type predict 
Koc values in the range of 230-240 MeV |42| . Moreover, 
it is possible to build bona fide Skyrme forces so that the 
SNM incompressibility is close to the relativistic value, 
namely 250-270 MeV. Hence, from the ISGMR experi¬ 
mental data, conclusion may be drawn that Ri 240 
± 20 MeV. Moreover, the constant of density dependence 
/3=1.5934±0.0085 fm^, which has the dimension of cross 
section for n = |, can be interpreted as the isospin av¬ 
eraged effective nucleon-nucleon interaction cross section 
in ground state symmetric nuclear medium. For a nu¬ 
cleon in ground state nuclear matter kp Ri 1.3 fm“^ and 
go ~ hkpc Ri 260 MeV and the present result for the ‘in 
medium’ effective cross section is reasonably close to the 
value obtained from a rigorous Dirac-Brueckner-Hartree- 
Fock [i^ calculations corresponding to such kp and go 
values which is Ri 12 mb. Using the value of the constant 
of density dependence /3=1.5934±0.0085 fm^ correspond¬ 
ing to the value of the parameter n = | along with the 
baryonic density of 0.1533 fm“^, the value obtained for 
the nuclear mean free path A is about 4 fm which is in ex- 
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cellent agreement [ij with that obtained using another 
method. 


VII. SUMMARY AND CONCLUSION 

In summary, the DDM3Y effective interaction which 
provides a unified description of elastic and inelastic scat¬ 
tering, proton-, a-, cluster- radioactivities and nuclear 
matter properties, also provides an excellent description 
of the /3-equilibrated neutron star matter which is stiff 
enough at high densities to reconcile with the recent ob¬ 


servations of the massive compact stars (25l - [^ while the 
corresponding symmetry energy is supersoft as preferred 
by the FOPI/GSI experimental data. The neutron star 
core-crust transition density, pressure and proton frac¬ 
tion determined from the thermodynamic stability condi¬ 
tion to be pt = 0.0938 fm“^, Pj = 0.5006 MeV fm“^ and 
Xp(i) — 0.0308, respectively, along with observed mini¬ 
mum crustal fraction of the total moment of inertia of 
the Vela pulsar provide a new limit for its radius. It is 
somewhat different from the other estimates [J and 
imposes a new constraint R > 4.10-I-O.SGM/Mq kms for 
the mass-radius relation of Vela pulsar like neutron stars. 
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